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We investigate chiral superconductivity which occurs in the electronic nematic state. A vortex 
state in a c-axis magnetic field is studied on the basis of the two-component Ginzburg-Landau 
model for nematic-chiral superconductors. It is shown that various vortex lattice structures are 
stabilized by nontrivial cooperation of nematicity and chirality in superconductors. In particular, 
the vortex lattice structural transition occurs when a square anisotropy parameter v is positive 
(negative) and the nematicity is induced along the [110] axis ([100] axis). We discuss nematic-chiral 
superconductivity in URu 2 Si 2 , Si^RuCH, and UPt 3 . An experimental test for the examination of 
nematic order and chiral superconductivity is proposed. 
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I. INTRODUCTION 

Recent studies on the strongly correlated electron sys¬ 
tems have explored the chiral superconductivity with 
broken time-reversal symmetry. A chiral ( p x ± ip y )-wave 
superconductivity analogous to the 3 He A phase has been 
established in S^RuO,* 1 ’ , and a chiral ( d zx ± id yz )-wave 
superconductivity in URu 2 Si 2 has been identified by re¬ 
cent experimental studies 3-7 . Spontaneous time-reversal 
symmetry breaking characteristic of chiral superconduc¬ 
tors has been observed in both compounds 5,6,8 ’ 9 . Fur¬ 
thermore, a recent polar Kerr rotation measurement de¬ 
tected broken time-reversal symmetry in UPt 3 10 , and 
thus a chiral superconducting (SC) state which belongs 
to the Ei u , E 2u , or E\ g representation is suggested. 

Interestingly, chiral superconductivity coexists with a 
nematic order in at least some of these superconductors. 
Nematic states in itinerant electron systems analogous to 
classical liquid crystals have been one of the highlights 
in recent condensed-matter physics 11 . A ’’quantum ne¬ 
matic liquid crystal” accompanying spontaneous rotation 
symmetry breaking has been studied in some strongly 
correlated electron systems. For instance, a nematic or¬ 
der arising from the fluctuating stripe order has been 
proposed for cuprate superconductors 12 , and a nematic 
state in a bilayer ruthenate Sr 3 Ru 2 07 has been investi¬ 
gated extensively 13 . Furthermore, a nematic fluctuation 
in Fe-based superconductors has been investigated 14,15 
and identified as a possible glue of Cooper pairs 16,17 . 

There is accumulating evidence for a nematic or¬ 
der in the so-called hidden-ordered state of URu 2 Si 2 18 . 
The magnetic torque 19 , cyclotron resonance 20 , NMR 21 , 
and x-ray scattering 22 measurements uncovered a bro¬ 
ken fourfold rotation symmetry below the hidden-order 
temperature T < ThOj although another NMR measure¬ 
ment did not detect any broken symmetry 23 . Since the 
SC transition is a continuous second-order phase transi¬ 
tion, the nematic order has to coexist with the supercon¬ 
ductivity below T c . Indeed, a signature of nematic order 
has been observed in the SC state 24 . A kink in the lower 
critical field H c i has been attributed to a second SC tran¬ 


sition due to the broken fourfold rotation symmetry. It 
is also known from early experimental results 25,26 that 
the sixfold rotation symmetry in the hexagonal crystal 
lattice of UPt 3 is broken in the SC state. The splitting 
of two transition temperatures and the existence of the 
A phase at zero magnetic field has been attributed to 
the effect of a nematicity 27,28 . Although the origin of 
broken rotation symmetry is still unclear, the antiferro¬ 
magnetic order 26,29 may cause the nematicity. What is 
clear is that a nematicity plays an essential role in the 
multiple SC phases in UPt 3 . Although a spontaneous 
nematic order does not occur in S^RuCH, a nematic¬ 
ity is artificially induced by uniaxial pressure. Indeed, 
an enhancement of the transition temperature due to the 
uniaxial pressure has been observed in S^RuCH 30 . Thus, 
all the chiral superconductors known up to now coexist 
with the nematicity. 

It is interesting to study the vortex state of these 
nematic-chiral superconductors because of the following 
two reasons. First, an effect of nematic order competes 
with the gradient mixing in order parameters due to a 
magnetic field. Although the nematicity favors a nonchi- 
ral state, such as the ( p x ±p y )-wave state or (d zx ± d yz )~ 
wave state, the magnetic field along the c axis stabi¬ 
lizes the chiral state, such as the (p x ± ip y )-wave state 
or (d zx ± id yz )-wave state, through the gradient mixing. 
Therefore, the order parameter in nematic-chiral super¬ 
conductors is nontrivial. Second, various vortex lattice 
structures are stabilized in the chiral superconductors be¬ 
cause of the gradient coupling of two-component order 
parameters 31-33 . Thus, it is expected that rich vortex 
lattice phases appear in nematic-chiral superconductors. 
From the other perspective, the vortex state in chiral 
superconductors is sensitive to the nematicity in the un¬ 
derlying electronic state. Therefore, the vortex lattice 
structure will be a sensitive probe for detecting the ne¬ 
matic order. 

In this paper we investigate the vortex state of 
nematic-chiral superconductors on the basis of the 
Ginzburg-Landau (GL) theory. Signatures of the ne¬ 
matic order in the chiral SC state are clarified. The or- 
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ganization of this paper is as follows. In Sec. II, we intro¬ 
duce the GL model and describe a variational method by 
which we study the vortex state. In Sec. Ill, we show the 
phase diagram for the order parameter and vortex lattice 
structure. Experimental results on URii2Si2, Sr2Ru04, 
and UPt3 are discussed and future experimental studies 
are proposed in Sec. IV. 

II. GINZBURG-LANDAU THEORY 

A. Ginzburg-Landau model 

In this section, we construct a two-component 
Ginzburg-Landau (GL) model for chiral superconduc¬ 
tors in the nematic state. Essential variables are 
two-component order parameters, (?7 a (r), r]b(r)), by 
which the pairing function is described as A(r, k) = 
’ 7 a(^) 0 a(fc) +r]b(r)(j) h (k). Here, </ a (fc) and 4 >b(k) stand 
for pairing functions in the momentum space which be¬ 
long to a two-dimensional representation of the crystal 
point group 34 . For instance, the chiral (p x ± ip y )-wave 
SC state in Sr 2 Ru04 belongs to the E u representation of 
the Z?4h point group 35 , while the chiral (d zx ± id yz )-wave 
state in URu2Si2 belongs to the E g representation. The 
chiral /-wave superconductivity, which belongs to the 
E2u representation of the D §h point group, has been sug¬ 
gested for a representative multicomponent superconduc¬ 
tor UPt3 27 . On the other hand, a recent thermal conduc¬ 
tivity measurement indicated the E\ n representation 36 . 
All of these chiral SC states are represented by two- 
component order parameters (?7 a (r), ?yb(r)) as above, al¬ 
though the pairing functions </> a (fc) and cf>b(k) depend on 
materials. 

First, we discuss a GL model in the absence of the 
nematicity. Assuming a weak nematicity, we later add a 
quadratic symmetry-breaking term. The symmetric part 
of the GL free-energy density has been obtained as 

F 0 = a(|?7a| 2 + |?7b| 2 ) + -^-(l?7a| 2 + I ^7 b | 2 ) 2 

+ y(f?af?b - c - c -) 2 + /?3|?7a| 2 |?7b| 2 

+ K l(l^x ? 7a|“ + |-Dy?7b| 2 ) 

T^d^b^ + I^al 2 ) 

+ K 3 [(D' x r] a )(D' y 1J h )* + C.C.] 

+ K 4 [(T> x r;b)(T , y ?7a)* + C.C.] (1) 

for chiral superconductors in the type-II limit 34 . Follow¬ 
ing the conventional notation, we denote a = ao(T /T 6 — 
1), with T 6 being the transition temperature at zero 
magnetic field, and covariant derivatives D/ = — i dj + 
( 27 t /<& 0 )Aj, with <P 0 = he/ 2 |e|. We omit gradient terms 
containing 2 derivatives because we focus on the vortex 
state in magnetic fields applied along the c axis. 

In the weak-coupling BCS theory, parameters are given 
by the following relations: fo/Pi = { 4 > 14 > 1 )ts / (^Dfs- 


/V/3i = 3 fo/Pi - 1 , k 2 /«i = (</> a u 2 ) Fg / (/> 2 u 2 ) FS , and 
K3/K1 = K 4/«1 = (4>a.^bV x Vy) FS / (<j)lvl) FS , where v x 
and v y are Fermi velocities in the ab plane and brackets 
(• • -) Fg denote an average over the Fermi surface. De¬ 
spite three independent parameters, P2/P1, k 2 /ki, and 
K3/K1, a single-parameter description of the GL model 
has been adopted in studies of chiral (p x ± ip y )-wave 
superconductivity 31-33 . As shown by Agterberg 31,32 , the 
coupling constants are represented by a single-parameter 
v as 37 

/VA = «2/ki = K3/K1 = (1 + v)/{Z - v) (2) 

when we assume pairing functions (/> a (fc), (j>\,(k)) = 
(v x (k), v y (fe)) (Agterberg model). Adopting simi¬ 
lar pairing functions for the chiral (d zx ± id yz )-wave 
SC state, (/ a (fe), </>b(fe)) = {v x (k)f(k z ), v y (k)f(k z )), 
with f(k z ) being an odd function (extended Agterberg 
model), we approximately obtain parameters shown in 
Eq. ( 2 ). Strictly speaking, the single-parameter descrip¬ 
tion breaks down since / 3 2 //?i 7^ k 2 /« 1 = K3/K1. How¬ 
ever, the deviation is negligible, P2/P1 — k 2 /ki 1, for 
a smooth function f(k z ). Thus, the single-parameter de¬ 
scription is also applicable to the chiral (d zx ± id yz )-wave 
superconductors. In crystals satisfying the D4 h point 
group, the parameter v indicates the square anisotropy 
in the Fermi surface, because v = 0 for the cylindrical or 
spherical Fermi surface, while v ^ 0 otherwise. On the 
other hand, the chiral F lu state in the Dgh point-group 
symmetry is similarly described as (</> a (fc), (f>b(k)) = 
(v x {k)f(k z ), v y {k)f{k z )) with f(k z ) being an even func¬ 
tion. Then, we obtain P2/P1 — k 2 /ki = K3/K1 = 1 / 3 , 
and thus, u = 0 in Eq. ( 2 ) irrespective of the Fermi 
surface 38 . Thus, we can rely on the single-parameter de¬ 
scription of the GL model for these chiral p- wave, d- wave, 
and /-wave superconductivities when the Fermi velocity 
and pairing functions have a smooth momentum depen¬ 
dence. Adopting the single-parameter description, we in¬ 
vestigate the vortex state in chiral superconductors. On 
the other hand, Eq. ( 2 ) seriously breaks down in the chi¬ 
ral E2u state in the Dq h point-group symmetry. Thus, 
the chiral F 2u state is beyond the scope of our study, but 
we will briefly discuss it in Sec. IV. 



FIG. 1. Primitive vectors of the vortex lattice ai and 0 . 2 . 
Two-dimensional coordinates ( x , y) and (x, y) and angles #l 
and 9 are illustrated. 

We rewrite the GL model using order parameters in 
the chirality basis, 771 = ^(77 a — irjb)/V^ and 772 = 7(?? a + 
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ir/b)/ a /2 with, 7 = (/?i/a 0 ) 1/2 and /?i = /3i - /? 2 + /V 2 = 
2/?i/(3 — ^). With the two-dimensional coordinate ro¬ 
tated through an angle 0 around the 2 axis (see Fig. 1), 
the dimensionless free-energy density is given by 

/° = ^rMM 2 + M 2 ) + \{\m\ A + M 4 ) 

+ 2|7i| 2 |7 72 | 2 -^[(r ?1 ^) 2 +c.c.] 

+ l-Dx^ll 2 + l-Dy^ll 2 + |-Dx%| 2 + | DyT]2 | 2 

X [(Ap7i)(Ap 72)* - (£>y??i)(Hy?72)*] +c.c.} 

+ i{i(e 2ie + ^e- 2ie ) 

x [(Ae^)^)* + (DymKDzm)*] + c.c.}, (3) 

where the unit of energy, length, and magnetic field are 
al/fii, £ = [(Ki + K 2 )/2a 0 ] 1/2 , and ^o/27t£ 2 , respectively. 
The covariant derivatives are denoted as D-j = —i<9y + 
A-r. We focus on the magnetic field along the c axis, 
and choose the vector potential A = — yHe In the 
reasonable parameter range \v\ < 1, a chiral SC state 
with ( 771 , 772 ) oc (1,0) or ( 771 , 772 ) oc (0,1) is stable at zero 
magnetic field. 

A weak nematicity leading to the violation of four¬ 
fold or sixfold rotational symmetry is taken into account 
by adding the symmetry-breaking term in the quadratic 
form 

hh = g{mih + c.c.) = ff7 2 (?7a77b + C.c.) (4) 

to fo- Thus, the total free-energy density is given by 

/ = fo + / 2 h. (5) 

The coupling constant g represents a lifting of degen¬ 
eracy between two orbital pairing functions 4> a (k) and 
</>b(fc) due to the nematic order. As a consequence of the 
symmetry-breaking term, double SC transitions occur at 
zero magnetic field. In the high-temperature phase near 
the SC transition temperature, T c2 < T < T c , a nonchi- 
ral state where ( 771 , 772 ) oc (1, ±i) is stabilized, while a chi¬ 
ral state is stabilized below T c2 . The time-reversal sym¬ 
metry is spontaneously broken in the low-temperature 
phase. The double SC transition has been observed in 
UPt 3 25 “ 28 , and a signature of double SC transition has 
been reported in URu 2 Si 2 24 . 

In Eq. (4) we assume a nematicity along the [110] direc¬ 
tion to be consistent with experimental observations in 
URu 2 Si 2 19 ~ 22 . The nematicity along the [100] axis is also 
investigated on the basis of our model. When we change 
the coordinate 9 —>• 9 -(- 7r/4 and take the phase factor 
( 771 , 772 ) —> (? 7 ie _17r / 4 ,772e 17r/,4 ) ) the fo term is almost in¬ 
variant except for the sign reversal v —> —v. Then, the 
symmetry-breaking term changes to 

S'i h = 9 (771772 + c.c.) = c/7 2 (|? 7 a 1 2 - | 77 b| 2 ) , (6) 


which is nothing but the symmetry-breaking term due 
to the nematicity along the [100] axis. A symmetry- 
breaking term of this form has been adopted for studies 
on multiple SC phases in UPt 3 27,28 . We change the sign 
of v instead of considering the symmetry-breaking term 
in Eq. (6) for studies of nematicity along the [100] axis. 

B. Variational method 

We investigate the order parameters and vortex lat¬ 
tice structure using the variational method. We assume 
variational wave functions of Cooper pairs so that the so¬ 
lution of the linearized GL equation is reproduced. First, 
we solve the linearized GL equation using the Landau- 
level expansion. Differentiating the quadratic terms in 
the fo term [Eq. (3)] with respect to 771 and 772 , we ob¬ 
tain the linearized GL equation in the absence of the 
symmetry-breaking term, 

/ 771 \ J_ / 1 + 2n+n_ e- 2ie U 2 _ - ve 2ie II 2 \ 

A V V2 )~ 11 V e 2W u% - ue~ 2id n 2 _ 1 + 2 n+n_ ) 

*(£)■ <7) 

where l c = 1/y/H and II± = —(l c /y/2)(D^ ± iDy). The 
solution for the minimum eigenvalue A = A m i n is repre¬ 
sented as 

( ipi+( r ) \ ( a in (e)ip 4n (r,p,a) \ f ) 

\lp 2 +(r) ) \a4n+2(9)ip4n+2(r,p,a) ) ’ ’ 

n>0 

where < p n (r,p,a) denotes the nth Landau-level wave 
function. The leading term is the lowest Landau 
level of the positive chirality component, and thus, 
(V>i+(7~),7/>2+(r)) ~ (a o (0)<po(r,p,a),0). We obtain an¬ 
other solution for the pairing state with a dominantly 
negative chirality, 

Vh -(r) \ _ ST ( b4n+2{9)(p4n+2(r,p,a) \ , . 

fp 2 -(r) ) hn{9)ip4n(r, p,a) )’ ^ ’ 

where the leading term is the lowest Landau level of 
the negative chirality component, (^i_(r), 't/> 2 — ( r ')) — 
(0,bo(9)tpo(r, p,a)). Coefficients a n {9) and b n (9) are nu¬ 
merically determined. We assume variational wave func¬ 
tions consisting of a linear combination of the two solu¬ 
tions: 



where |(7+| (|C_|) represents the weight of Cooper pairs 
having dominantly positive (negative) chirality. The vari¬ 
ational wave function is justified near the transition tem¬ 
perature and for a small symmetry-breaking term g, al¬ 
though the reconstruction of higher Landau levels affects 
the vortex lattice structure at low temperatures 33 . Our 
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main result is concerned with the vortex lattice struc¬ 
tural transition near T c , and thus, the variational wave 
function is appropriate. 

In order to study the vortex lattice structure, we adopt 
a general form of the nth Landau-level wave functions 39 


<Pn(r,p,c r) 


V2 n 7r 1 / 2 n! 


E 


Cm C 


27ri(m—l/2)x/a 


x H n 


^ y__yrn}j e -{y-y m ) 2 /2ll ^ 


( 11 ) 


where c m = e 17rm ( p+1_rrap \ y m = l c \j2^ii(m — 1/2), and 
H n (y) is the Hermit polynomials. The vortex lattice 
structure is determined by the variables p = (b/a) cos$l 
and <j = (b/a) sin^L- Primitive vectors are a 4 = ae^ and 
a 2 = 6(cos+ sin0 L ey) = a(pe x + erey) (see Fig. 1). 
The area of the unit cell is \ai x a 2 \ = a&sinfA = 27 tZ 2 . 
The rectangular and centered rectangular lattices are 
formed for p = 0 and p = 1/2, respectively. The 
square and triangular lattices are special cases of them; 
(p, cr) = (0, 1) or (p, cr) = (1/2,1/2) in the square lattice, 
and (p, cr) = (1/2, a/3/2) in the triangular lattice. The 
angle of a primitive vector a 4 from the original x axis is 
0, and it is dealt with as a variational parameter. Thus, 
variational parameters for the vortex lattice structure are 
p, <7, and 6. 

Substituting Eqs. (8)-(ll) into Eqs. (3) and (4) and 
integrating over the unit cell, we obtain the free-energy 
density as 


F(C+, C-,p, a, 9) = </> uc = ^ j f(r)d 2 r. (12) 

The brackets (• • -) uc denote an average over the unit 
cell. Optimizing variational parameters (C+, C-, p, a, 6) 
to minimize the free-energy density, we determine the 
order parameters and vortex lattice structure. 


III. CHIRAL SUPERCONDUCTING STATES 

A. Vortex lattice in the non-nematic state 

In this section, we investigate the vortex state in the 
absence of nematicity. Thus, we consider a tetragonal or 
hexagonal system and choose the parameter g = 0. The 
effects of nematic order on the chiral SC state are studied 
in the following sections. 

First, we determine the vortex lattice structure near 
the upper critical field. Since the symmetry-breaking 
term / 2 h is absent, a solution of the linearized GL equa¬ 
tion ( 771 , 772 ) = C+( 7 / 1 - 1 - (r), ip 2+ (r)) minimizes the free 
energy near the SC transition. Thus, C- = 0, and the 
vortex lattice structure is determined so as to minimize 
the Abrikosov parameter /? A = 2 (/ 4 ) uc /((| t/i + | 2 ) uc + 
(|t/ 2 -i-| 2 ) uc ) 2 , where / 4 is the quartic term in Eq. (3). Fig¬ 
ure 2 shows the Abrikosov parameter for various vortex 



FIG. 2. (Color online) Abrikosov parameter /3 a for the rect¬ 
angular lattice (p = 0, red solid line) and centered rectangular 
lattice (p = 1/2, black long-dashed line) as a function of v. 
The variational parameter a is optimized to minimize /3 a- The 
Abrikosov parameter for the square lattice [(p, a) = (0,1)] 
is shown by the green short-dashed line. The other varia¬ 
tional parameter 6 is fixed to 9 = 0 (9 = 7t/ 4) for v > 0 
(v < 0), for which the Abrikosov parameter is optimized 
when \v\ > 0.0114. When \u\ < 0.0114, the optimized angle 
is 6 = 71/4 (6 = 0) for v > 0 (v < 0), but the 6 dependence of 
/3 a is negligible. 
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FIG. 3. Schematic of the vortex lattice structure near the 
upper critical field. We show the result for v > 0. The vortex 
lattice is rotated 7r/4 for v < 0. 


lattice structures as a function of the anisotropy parame¬ 
ter v 1 and wc illustrate the vortex lattice structure which 
minimizes the free energy in Fig. 3. It is shown that the 
centered rectangular lattice, square lattice, rectangular 
lattice, and, again, centered rectangular lattice are sta¬ 
bilized with increasing \v\. These results are consistent 
with previous works which investigated the vortex lattice 
structure in the small- \u\ region 31-33 . Agterberg showed 
that the square lattice is stable unless the anisotropy pa¬ 
rameter is extraordinary small, \v\ < 0.0114 31 - 32 . The 
triangular lattice is formed at v = 0 like in conventional 
superconductors, and it is deformed with increasing \u\. 
The parameter a decreases from \/?>/2 to 1/2 when \v\ 
increases from 0 to 0.0114. The qualitatively same results 
have been obtained by Kita using a sophisticated calcu¬ 
lation which takes into account higher Landau levels and 
screening current 33 . He also found that the rectangu¬ 
lar lattice is stabilized for a large anisotropy parameter 
\v\. Our calculation reproduces these results and, fur¬ 
thermore, shows that the centered rectangular lattice and 












5 


rectangular lattice are stabilized for 0.27 < \v\ < 0.58 
and 0.58 < \v\ < 0.7, respectively. We do not consider 
a further large parameter \p\ > 0.7 since the numerical 
convergence becomes worse. 

Next, we show the phase diagram against magnetic 
fields and temperatures (H-T phase diagram) for two 
anisotropy parameters v = 0.15 and v = 0.3. For 
v = 0.15, the square vortex lattice is stable near the tran¬ 
sition temperature T = T C (H ), as illustrated in Fig. 3. 
Figure 4 shows that the structural transition occurs at 
a moderate temperature T = T C 2 (H), and the rectan¬ 
gular vortex lattice is stabilized below T c2 (F7). Figure 
5 shows the magnetic field dependences of variational 
parameters |C±| and the ’’nematicity of the vortex lat¬ 
tice” b/a at T/T c ° = 0.4. It is shown that the emer¬ 
gence of nematicity b/a — 1 coincides with CL, although 
C- = 0 and b/a = 1 in the high magnetic field region, 
H/H C 2(0) > 0.38. Thus, the structural transition in the 
vortex lattice accompanies the mixing of chirality in the 
order parameter. 



FIG. 4. (Color online) Phase diagram of the vortex lattice as 
a function of the dimensionless magnetic field H/H c 2 ( 0 ) and 
temperature T/T® for (u,g) = (0.15,0). The red solid line 
shows the SC transition temperature T C (H), while the black 
dashed line shows the second transition temperature T C 2 (H). 
The vortex lattice structures are shown schematically. The 
phase diagram is independent of the sign of u, while the vortex 
lattice for v < 0 is rotated 7 t/4 from that for v > 0. 


On the other hand, the vortex lattice structural transi¬ 
tion does not occur for v = 0.3. Then, the centered rect¬ 
angular lattice is formed in the whole SC state (Fig. 6 ). 
Because of the orthorhombic symmetry of the vortex lat¬ 
tice, the order parameter for pairing with dominantly 
negative chirality is finite, CL 7 ^ 0. However, the pos¬ 
itive chirality is favored by the linear coupling of mag¬ 
netic field and chirality in Cooper pairs 31 , and therefore, 
|C+| |C_|. A similar H-T phase diagram is obtained 

for the isotropic case v = 0. Then, the triangular vortex 
lattice is stabilized. 




FIG. 5. (Color online) Magnetic field dependence of (a) order 
parameters |C±| and (b) a structural parameter b/a (= a) 
at T/T c ° = 0.4 for (v,g) = (0.15,0). The other structural 
parameters are #l = 7 t /2 (p = 0) and 6 = 0. 



FIG. 6. (Color online) Phase diagram of the vortex lattice 
for (v,g) = (0.3,0). As shown schematically, the centered 
rectangular lattice is stabilized in the whole SC state. 


B. Vortex lattice in the nematic state (u > 0) 

Now we turn to the main topic of this paper. We study 
the vortex state of chiral superconductors which coexist 
with the nematic order. The GL model with a finite 
symmetry-breaking term (g 7 ^ 0) is analyzed. The sign 
of the coupling constant g is not important at all, and 
thus, we assume g > 0. For g < 0, the vortex lattice 
rotates tt/2. We choose g = 0.05 unless we explicitly 
state otherwise. The SC double transition occurs at zero 
magnetic field, and the splitting of two transition tem¬ 
peratures is 10% of T c ° for our choice of the coupling 
constant g. This splitting is consistent with the experi¬ 
mental data indicating a double transition in URu 2 Si 2 24 
and UPt 3 25-28 . 
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In a magnetic field H > 0, the chiral SC state with 
positive chirality ( 771 , 772 ) oc (1,0) is favored by the gra¬ 
dient mixing in order parameters 31 . Then, the second 
SC transition is smeared and it changes to the chiral- 
nonchiral crossover (C-NC crossover in Figs. 7, 11, and 
12). We show here that the vortex lattice structural tran¬ 
sition occurs as a result of the C-NC crossover when the 
anisotropy parameter is positive, v > 0. Owing to the 
symmetry-breaking term, the phase diagram is no longer 
independent of the sign of v. Thus, we study the case 
with v > 0 in this section, and the case with v < 0 is in¬ 
vestigated in the next section. The isotropic case, v = 0, 
will be discussed for UPt 3 in Sec. IV. 


i, r 11 hi 



FIG. 7. (Color online) Phase diagram for {v, g) = (0.15, 0.05). 
The red solid line shows the SC transition temperature T C (H). 
The green squares and blue triangles show the critical point 
of the vortex lattice structural transition. The black dot¬ 
ted line depicts the C-NC crossover, which we define by 
|C_| = 0.7|C+|. Phase I’ is the nonchiral SC state in which 
the pairing function is approximately described as A(fe) ~ 
4>a,(k) — 0b(fc), while phases I, II, and III are the chiral SC 
states in which A(fc) ~ </> a (fc) + i(/>b(fc). The vortex lattice 
structures in phases I (I’), II, and III are illustrated on top of 
the phase diagram. 


First, we show the phase diagram for v = 0.15 in Fig. 7. 
Indeed, the vortex lattice structural transition occurs, 
and three phases, I (I’), II, and III, appear. Although 
the square or rectangular lattice is formed in the absence 
of the nematicity (see Fig. 4), a centered rectangular lat¬ 
tice is stabilized in phases I and T near T c . This is intu¬ 
itively understood as follows. Although a finite chirality 
is slightly induced by the gradient mixing, the order pa¬ 
rameter is approximately nonchiral, ( 771 , 772 ) oc (1,— i), 
in the high-tenrperature region near T c . Thus, an elon¬ 
gated triangular vortex lattice is stabilized as in single¬ 
component superconductors. Note that the elongated tri¬ 
angular lattice is equivalent to the centered rectangular 


lattice. 

We would like to stress that the symmetry-breaking 
term significantly affects the vortex lattice structure even 
below the C-NC crossover temperature. We define the 
C-NC crossover line dividing phases I and I’ using of the 
variational parameters as |C_| = 0.7|C+|. We see that 
the centered rectangular lattice (which equals the elon¬ 
gated triangular lattice) is stabilized well below the C-NC 
crossover temperature. Thus, the vortex lattice structure 
in chiral superconductors is sensitive to the nematic or¬ 
der. This is because of the small stiffness of Abrikosov 
vortex lattice and the large coupling of nematicity and 
chirality in superconductors. 

The vortex lattice structure in phase II is interpreted to 
be a distorted rectangular lattice. A nematicity gives rise 
to a stress along the [ 110 ] axis, and therefore, the rectan¬ 
gular lattice is deformed to an oblique lattice (p 7 ^ 0 , 1 / 2 ) 
in phase II. In particular, a marked effect of nematicity 
appears near the second-order phase transition to phase I, 
as we show the temperature dependence of variational pa¬ 
rameters in Fig. 8 . The structural parameter is p = 1/2 
in phase I (centered rectangular lattice) and decreases to 
p = 0 (rectangular lattice) with decreasing temperature 
in phase II. At the same time the angle of the primitive 
vector 9 rotates from 0 = 7 r/ 4 to 6 | ~O. The second- 
order transition from phase I to phase II is accompanied 
by the spontaneous violation of the reflection symmetry 
along the [110] axis. Thus, we see marked signatures 
of nematic order in the vortex lattice structure in the 
low-temperature phase II. We show typical vortex lat¬ 
tice structures in phases I and II in Fig. 9(a) and 9(b), 
respectively. 

The square vortex lattice is deformed to the centered 
rectangular lattice in the high-field phase III owing to 
the symmetry-breaking term. Then, 9 l 7 ^ 7t/2; however, 
the deformation is negligible, as shown in Fig. 10. The 
structural phase transition between phase II and phase 
III is characterized by the increase in 6 /a, as it occurs in 
the absence of nematicity (see Fig. 5). 

Although the phase transition between phases II and 
III is specific for the parameter v = 0.15, the vortex 
lattice structural transition between phases I and II is 
ubiquitous as it is induced by the C-NC crossover. In¬ 
deed, the latter occurs for a wide range of the anisotropy 
parameter v. For instance, we show the phase diagram 
at v = 0.3 (Fig. 11). It is shown that the vortex lattice 
structural transition is induced by the C-NC crossover. 
The centered rectangular lattice is stabilized in the high- 
temperature phases I and I’, while the oblique lattice is 
stabilized in the low-temperature phase II, similar to the 
case of v = 0.15. Although the centered rectangular lat¬ 
tice is stable in the absence of the nematicity for v = 0.3 
(see Fig. 6 ), it is deformed to the oblique lattice because 
the primitive vectors a\ and 0,2 are not parallel to the 
stress along the [110] axis. The orientation of the vor¬ 
tex lattice is significantly rotated from 9 = 7t/ 4 to 9 ~ 0 
by decreasing the temperature in phase II. Interestingly, 
a structural parameter p shows a marked change from 
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FIG. 8. (Color online) Temperature dependence of (a) order parameters \C±\ and structural parameters (b) a, (c) p, and (d) 
6 at H/H c2 ( 0) = 0.1 for {u,g) = (0.15, 0.05). 


(a) \r]( r ) 



(b) 17 /(r) 



FIG. 9. (Color online) Vortex lattice structures (a) in phase 
I and (b) in phase II for (u, g) = (0.15,0.05). We choose 
T/T° = 0.8 and H/H c 2 ( 0 ) = 0.1 in (a) and T/T° = 0.5 and 
H/Hc 2 ( 0 ) = 0.1 in (b). We show the amplitude of the order 
parameter |p(r)|, which is obtained by \r/(r)\ 2 = |J 7 i(r)| 2 + 

Mr)| 2 - 


p = 0.5 to p ~ —0.4 at the same time. 

Our results for v = 0.15 and v = 0.3 imply that the 
angle of a primitive vector from the x axis 9 plays an 
important role in the vortex lattice structural transition. 
The x axis is no longer a principal axis of the electronic 
state for j ^ 0. Because the primitive vector should be 
parallel to the principal axis in the nonchiral state, we 
obtain 0 = ±7t/ 4 above the C-NC crossover temperature. 
The vortex lattice structural transition is induced by the 
C-NC crossover when we obtain 9 ~ 0 in the chiral state. 
As shown in Fig. 3, this condition is satisfied for v > 
0.00114 since the vortex lattice structure in the chiral 
state is little affected by the symmetry-breaking term. 
As expected from this consideration, the vortex lattice 
structural transition between phases I and II disappears 
when v < —0.00114. Then, 9 = ±7r/4 in the whole SC 
state, as we will show in the next section. 


C. Vortex lattice in the nematic state (v < 0) 

As expected from the above discussions, Fig. 12 shows 
the H-T phase diagram without any structural phase 
transition for v = —0.15 and v = —0.3. The C-NC 
crossover occurs, as indicated by the dotted lines. How¬ 
ever, the vortex lattice structural transition does not 
occur because the stress due to the nematicity is ap¬ 
plied parallel to a primitive vector of the vortex lattice 
(9 = 7t/ 4). The rectangular lattice and centered rect¬ 
angular lattice are stabilized in the whole SC state for 
v = —0.15 and for v — —0.3, respectively. The effect of 
nematicity is only the increase in the parameter a. The 
sign of the coupling constant g determines the orienta- 





















FIG. 10. (Color online) Magnetic field dependence of (a) order parameters \C± \ and structural parameters (b) b/a, (c) #l, and 
(d) 9 at T/T° = 0.1 for (v,g) = (0.15,0.05). We show b/a and 0 l instead of a and p. 




(a) (v, g)=(-0.15, 0.05) 


(b) (v, g)=(-0.3, 0.05) 


FIG. 11. (Color online) Phase diagram for ( v,g ) = (0.3, 0.05). 
The order parameter and vortex lattice structures in phases 
I, I’, and II are similar to those in Fig. 7. 


FIG. 12. (Color online) Phase diagram for v < 0. (a) The 
rectangular lattice (p = 0,9 = 7r/4) is stable in the whole 
region of the SC state for (v, g) = (—0.15, 0.05), while (b) the 
centered rectangular lattice is stable for (v,g) = (—0.3,0.05). 


tion of the elongated vortex lattice. When we choose 
a > 1 for p = 0 and er > 1/2 for p = 1/2, the orientation 
is 9 = 7t/ 4 for g > 0, while 9 = —7t/ 4 for g < 0. For 
v = —0.15, the square-rectangular structural transition 
occurs at g = 0, but it is smeared by the external nematic 
order. The square lattice is deformed to the rectangular 
lattice owing to nematic order. 


IV. SUMMARY AND DISCUSSION 

In this paper we studied chiral SC states coexisting 
with a nematic order. The phase diagram of the order 
parameter and vortex lattice structure has been clarified 
on the basis of the two-component GL model. It is shown 


that the vortex lattice structure is sensitive to the chi¬ 
rality in Cooper pairs as well as to the nematicity in the 
electronic state. Because the chirality and nematicity co¬ 
operate in a nontrivial way, various vortex lattice struc¬ 
tures are stabilized. In particular, the structural phase 
transition occurs in the vortex lattice when the nematic¬ 
ity is along the [110] axis and the anisotropy parameter 
is positive v > 0, or when the nematicity is along the 
[100] axis and v < 0 37 . Otherwise, the vortex lattice 
structural transition is not induced by the nematic or¬ 
der, although the structural parameters are affected by 
the C-NC crossover in the pairing function. The struc¬ 
tural phase transition and structural change discussed 
in this work are distinguished from those due to the 
anisotropy in the Fermi surface and SC gap 40,41 , gra- 
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dient mixing in several irreducible representations 42 , and 
the Pauli depairing effect 43,44 . The former occurs in the 
low-magnetic-field region near T c , although the latter ap¬ 
pears in the high-magnetic-field region. 

Finally, we discuss the vortex state in the chiral su¬ 
perconductors on the basis of our results. Our study 
has been mainly focused on a heavy-fermion supercon¬ 
ductor URri 2 Si 2 18 in which nematic order has attracted 
much attention recently. The nematicity observed by sev¬ 
eral experiments 19-22 is along the [110]-direction. When 
we assume that two heavy Fermi surfaces around the 
M point cause the superconductivity, the anisotropy pa¬ 
rameter is estimated to be v = 0.41 45 on the basis of 
the band structure calculation 46 and the Shubnikov-de 
Haas measurement 47 . Thus, it is reasonable to assume 
a positive v, although the anisotropy parameter should 
be affected by the multiband structure. According to our 
results, it is expected that the structural phase transition 
will occur in the vortex state under the c-axis magnetic 
field. If it were be observed by a small-angle neutron 
scattering (SANS) measurement 48 , for instance, clear ev¬ 
idence for both nematic order and chiral superconductiv¬ 
ity would be obtained. 

SANS measurements reported a square vortex lattice 
in Sr 2 Ru 04 49 . The orientation is 9 = 7 r/ 4 , implying a 
small negative parameter v < 0 , and therefore, a ne¬ 
maticity along the [ 100 ] axis may induce the vortex lat¬ 
tice structural transition. Then, the symmetry-breaking 
term deforms the square lattice to the centered rectan¬ 
gular lattice at low temperatures. On the other hand, 
the nonchiral SC state with another centered rectangu¬ 
lar vortex lattice is stabilized near T c . The two vortex 
lattices have the same symmetry; however, the structural 
parameters are quite different. For instance, our calcu¬ 
lation for v = —0.05 and g = 0.01 shows that the struc¬ 
tural parameter is 2er ~ 1.1 in the chiral SC state but 
2a ~ 3.0 in the nonchiral SC state. Thus, the vortex lat¬ 
tice structural change accompanying the C-NC crossover 
may be caused by the uniaxial pressure along the [ 100 ] 
axis. Interestingly, a significant enhancement of the tran¬ 
sition temperature due to the uniaxial pressure has been 
observed 30 . However, we should discuss the spin degree 
of freedom in Cooper pairs because Sr^RuCH is consid¬ 
ered to be a spin-triplet superconductor 1,2 . When the 
spin-orbit coupling in Cooper pairs is small enough to al¬ 
low the rotation of the d vector, a helical SC state, such 
as d = p x x — p y y , is stabilized in the c-axis magnetic 
field 50 . This state is beyond the scope of this paper, but 
it is expected that the nematicity affects the vortex lat¬ 
tice structure through the crossover in the pairing func¬ 
tion. For instance, the nematicity along the [100] axis 


stabilizes a planar state d = p x x or d = p y y near T c , 
while the helical state is robust at low temperatures. 

The symmetry of the pairing state in UPt 3 is 
controversial despite intensive studies for more than 
three decades. A recent polar Kerr rotation mea¬ 
surement found a spontaneous time-reversal symmetry¬ 
breaking, which implies a chiral SC state d(k) = 
[^a(fe) ± i 0 b(fc)[ c 10 , consistent with the suppression of 
the upper critical field along the c axis 2 '. On the 
other hand, NMR data indicate a rotation of the d 
vector 51 and implies a helical SC state, such as d(k) = 
</> a (k)b + <j>b(k)c. As for the orbital symmetry, the E 2u 
and Ei g representations have been considered as pre¬ 
dominant candidates 27,28 . On the other hand, a recent 
thermal conductivity measurement is consistent with the 
E\ u orbital symmetry with a hybrid nodal structure 36 . 
Among these SC states, the chiral Ei u state and chiral 
Ei g state are described by the extended Agterberg model, 
and then, the square anisotropy parameter has to be zero, 
v = 0 , owing to the hexagonal symmetry in the crystal 
lattice. Unfortunately, this model seems to be incompat¬ 
ible with several experimental results. First, the SC dou¬ 
ble transition is smeared in the c-axis magnetic field, al¬ 
though it has been observed in experiments 25-28 . Second, 
the vortex lattice structure significantly changes through 
the C-NC crossover, although only a realignment without 
any structural deformation has been observed by a SANS 
measurement 52 . In order to avoid these discrepancies, we 
have to assume parameters in the GL model without re¬ 
lying on the extended Agterberg model. The parameters 
allowed by symmetry 38 , k 2 — «i and K 3 = K 4 ~ 0, may 
be compatible with those experiments and may be con¬ 
sistent with the nearly isotropic upper critical field in the 
basal plane 36 . Although these parameters are naturally 
obtained by the weak-coupling theory for the E 2u state 27 , 
a fine-tuning of pairing functions and Fermi surfaces is 
required for the chiral Ei u and Ei g states. 
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